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Abstract 

The magnetic analog of odd-time superconducting order is introduced. It 
is shown that the classification of possible odd-time magnetic states admits 
the existence of a novel state - the chiral spin nematic - which is character- 
ized by an odd-in-time spin-spin correlation function. The known example 
of the chiral spin liquid appears naturally in this approach. The equal-time 
correlation function of these odd-in-time states involves non-trivial three spin 
correlations. 
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The question of broken symmetry states characterized by some correlation function which 
is odd in time, was revived recently in the context of superconductivity. An odd-time, 
or equivalently, odd-frequency, superconducting ordering was introduced in the pioneering 
work of Berezinskii who considered spin-triplet odd-time pairing as an alternative to 
the conventional theory of 3 He pairing. Interest in odd-time ordering was renewed after the 
work of Balatsky and Abrahams ||, in which the extension of Ref. [1] to odd-time pairing 
for spin-singlet superconductors was given. 

In this paper, in analogy to odd-time superconductivity, we introduce the concept of 
magnetic states with odd-time correlations. We shall discuss the symmetry constraints for 
such states and shall classify the different possibilities. 

For odd-time superconductivity, the possibility of odd-time states follows immediately 
from the basic symmetry constraint F(k, t) = F(— k, —t) for the general spin-singlet 
anomalous correlator F(k,t) = i((J y ) a p( T t Ck a (t) c_k/3(0)) . This symmetry constraint is an 
immediate consequence of the Fermi statistics of the Ck« operators and of the singlet nature 
of the anomalous correlator. The consideration for spin-triplet pairing is analogous and was 
done by Berezinskii |T|. Then one finds two classes of states: 1) even-time pairing (BCS) 
states and 2) odd-time pairing states with F(k, — t) = — F(k, t). The orbital parities of these 
states are opposite [even(odd) and odd(even) for singlet (triplet) respectively] as follows from 
the symmetry equation. 

The generalization of odd-time ordering to spin systems is straightforward and requires 
the symmetry equation for the dynamic correlation function of the spin density <Si(r, t) where 
i = 1,2,3 denotes cartesian components. The spin-spin correlation function is Ay(r, r'\t) = 
(T t Si(r,t)Sj(v',0)). The symmetry equation which follows directly from the properties of 
the time-ordering operator in the definition of Ay is A^-(r, r'|t) = Ajj(r',r| — t). In matrix 
form, this is: 

PT A(r,r'|t)=A(r,r'|t), (1) 

where P is the spatial parity operator r — r' — > r' — r, T is the time inversion and A is the 
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transposed matrix Ay = A^. This equation, valid for any rank spin S, allows the existence 
of odd-time magnetic states. These states are characterized by a nontrivial time structure 
of the spin-spin correlation function with Ay(r, r'\t) being an odd function of time. 

The main purpose of this paper is to classify magnetic states with odd-time magnetic 
correlations. We shall give the symmetry analysis for possible odd-time magnets. It follows 
closely that for even-time magnets which have nontrivial spin correlation functions. All the 
known chiral states are recovered in this approach: The chiral spin liquid (CSL) and the 
ferromagnet (a rather dull example of a chiral state) are among them. A new magnetic state 
is found as well: It is a state which is the odd-time analog of the spin nematic. The spin 
nematic was first considered by Andreev and Grishchuk || . The new odd-time state is also 
characterized by nematic ordering in spin space ||, but with broken T and P. We call this 
state a "chiral spin nematic" (CSN). We shall discuss the physical properties of this state 
and we shall show that the low energy lagrangian is identical to that for a conventional spin 
nematic [[|. 

We begin with the symmetry of the magnetic correlator. We consider states in which 
spontaneous breakdown of the 0(3) spin rotation group occurs without the appearance of 
an average microscopic spin density so that (Si(r,t)) = 0. This restriction leads us to 
consider a spin-spin correlation function Ay(r, r'\t), which is time-dependent and in general 
is a function of two coordinates r, r'. As shown by Andreev and Grishchuk ||, spin-spin 
correlations without time-reversal violation can still describe a nontrivial magnetic order, 
i.e. a spin nematic. In that case, Ay(r,r'|0) transforms as a tensor representation of 0(3). 
This will also be the case for odd-time correlations described by Ay for t ^ 0. 

Following the analysis of the Andreev and Grishchuk ||, we consider the invariance 
group of the system. It is given by the product 0(3) x T x G, where 0(3) is the spin 
rotation group, G is the space group of the lattice and T is time reversal, equivalent to spin 
inversion. G defines the symmetry of spin scalars; it contains the crystal point group and 
the translations of the lattice. 
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Aij is a 3 x 3 matrix containing 9 elements. It can be decomposed as a sum of terms 
which transform as particular representations of 0(3). From the decomposition 9=1+3+5 
we find that I = 0, 1, 2 can be present. Therefore we write: 

A^r, r'\t) = A(r, r'|t)<% + e ijk B k (r, r'\t) + Q -(r, r'\t). (2) 

where A(r,r'\t) is a scalar (7 = 0) and is the only term which contributes to the trace of 
Ay(r, r'\t). Bk(r,r'\t) is a vector (/ = 1) and Qij(r,r'\t) is the tensor part (/ = 2). The 
latter is symmetric and traceless. The quantities A, B k and Qij describe a spin liquid 
with no ordering, magnetic ordering with a pseudovector order parameter and spin nematic 
ordering, respectively. This is the case whatever the parity and time reversal properties of 
these objects. 

For the case of the even-time states, the order parameter can be found from the equal- 
time correlator Ajj(r,r'|t = 0). Then we recover the classification of Ref. ||. In the even- 
time case, Ajj(r,r'|t) in Eq. (|2|) is even under time reversal. Then Bk(r,r'\t) and Qij(r,r ; \t) 
describe two different nematic states, the so-called p and n nematics respectively [Q]. The 
A-term describes a spin liquid with no spin ordering. Complete details for the even-time 
case may be found in Ref. 0. 

We now turn to odd-time correlations. In the case of odd-time superconductivity, it is 
well-established that the equal-time correlators which describe the condensate contain 
composite operators which involve a conventional Cooper pair bound with, for example, a 
spin-density operator. In the magnetic case, we shall implement exactly the same approach 
and the odd-time magnetic correlators will describe the condensate of the higher order spin 
operators. For any Ay(r,r'|£) which is odd in time, the time derivative will be nonzero 
(we assume analyticity of Ay at small t) at t — and represents a thermodynamic order 
parameter for a composite condensate: 

d t K 3 {v,v'\t) t= , = (T t [d t 3 i (r,t)]3 j {T f ,0))^ (3) 

The equation of motion for Si can be used to find the general result 
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d t Si(T,t)=i[H,S i (T,t)] = e Um Si(T,t)M m (T,t), (4) 

where M m (r, t) is the "molecular field" for the hamiltonian TC. If the hamiltonian is bilinear 
in spin operators, the general form of the molecular field will be 

M m (r) = J dr'K mn (r,r')S n (r'). (5) 

The coupling K mn (r,r') is an explicit function of the two co-ordinates r, r'. Finally, we 
obtain 

flfcAyfo r'|*) t=0 = e llm J dr"K rnn (v, ^(S^S^S^r')) . (6) 
If the hamiltonian has the form 

U = - £ / drdr'S m (r)L mn {r, r')S n (r'), (7) 

mn 

then the molecular field kernel is given by 

K mn (r, r') = L mn (r, r') + L nm (r', r) = 2L mn (r, r'). (8) 

We see that any pure pairwise exchange terms of the hamiltonian will not contribute to the 
time derivative of odd-time correlators, since the corresponding kernel K mn (r, r") for pairwise 
exchange is symmetric in r, r", while the correlator on the r.h.s of Eq. (f|) is antisymmetric. 
Therefore for a non-zero result, K(r, r") must contain a spatially-odd component. We shall 
give an explicit example of a hamiltonian which produces these equations of motion for the 
specific case of a CSL. In general, complicated interactions are required for chiral magnetic 
order. 

We now discuss the specific phases which can be found in odd-time magnets. The possible 
terms in the decomposition Eq. (Q) will be considered in turn. 

1). Chiral Spin Liquid. Consider a spin-spin correlator Ay containing only the A-term 
in the r.h.s. of Eq. (0). The spin scalar term A(r,r'\t) is odd in time and, therefore is 
odd under spatial parity P, as can be easily seen from the symmetry equation ([I]) for any 
cartesian-symmetric odd-time A^-. Taking the time derivative as in Eqs. (3-6), we find 
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d t A(r,r'\0) = J dr"K(r,r")(S(r) x S(r") ■ S(r')>, (9) 

where, for the assumed CSL state with no anisotropy in spin space, we have taken 
K mn (r,r") = 5 mn K(r,r"). Note that d t A(r, r'|0) obtained in this way is explicitly odd 
under P. 

The necessary condition for the odd-time CSL phase is the existence of a non-zero corre- 
lator Xp(rir 2 r 3 ) = (S(ri) x S(r 2 ) • S(r 3 )), which when defined on a plaquette is exactly the 
usual CSL order parameter |§. We argue, therefore, that an alternative way to describe a 
CSL phase is in terms of an odd-time spin singlet correlator A(r, r'\t) = (T t Si(r, t)Si(r', 0)). 
The well-known picture of the CSL state as the state where all spins precess so that local 
handedness is preserved is naturally supported in this description. In order for the time 
derivative dtA(t = 0) to exist, the coupling K(r, r') should contain a part which is odd in 
r — r'. This puts a severe constraint on the possible exchange models which can support the 



CSL ground state. One of them is given by the hamiltonian [10 



« = "4 E [Si x S 2 • S 3 ] Pi [S 4 x S 5 • S 6 ] P . , (10) 

with A-positive. The sum runs over nearest neighbor plaquettes Pi, Pj with spins 1,2,3 
and 4, 5, 6 belonging to the plaquettes Pj and Pj respectively || and we consider here the 
triangular lattice. This hamiltonian has a CSL ground state with ferromagnetic ordering 
of chiralities with nonzero Xp(ri, r 2 , r 3 ). It can be shown that the equations of motion 
obtained from this hamiltonian do indeed lead to Eq. (||) for d t A, where the spin scalar 
kernel K(r,r') is a certain combination of the pair spin (S(rj)S(rj)) correlation functions. 

2). Ferromagnet. Consider the Bk term in Eq. (0). B>k transforms as a vector under the 
spin rotation group 0(3); it is even under spatial inversion P and odd under time reversal T, 
as follows from Eq. ([I]). A pairwise exchange interaction hamiltonian, therefore, is sufficient 
for the present discussion. One can take TC — | / drdr'K(r,r') S(r) ■ S(r'). Following the 
same approach as for the CSL case, we find for the time derivative d t B> 

<9 t 5(r,r'|0) = J dr"K(r,v"){([S(r") ■ S(r')]S(r)> - S 2 (S(r")>}. (11) 



Therefore, we conclude that a nonzero P-term in the decomposition of Eq. (^) describes 
a vector ordering (S(r)) 7^ 0. This is trivial ferromagnetic order. In particular, for r = 
r' it follows from Eq. ( p]) that dtB(r, r|0) oc (S(r)). The ferromagnet corresponds to 
the reducible three-spin correlator in Eq. ([]) and can be described in terms of a single 
spin operator expectation value. This contradicts our initial restriction to states having no 
average spin density and we do not consider this ferromagnetic case any further. 

3). Chiral Spin Nematic. As remarked earlier, the last term in the decomposition of 
Eq. (pj) represents a spin nematic order belonging to the tensor (7 = 2) representation 
of 0(3). The crucial difference from previously considered spin nematics [0 is that in 
our case Qij(r, r'\t) is an odd function of t and therefore it is an odd function under P. 
In the presence of an external magnetic field, which expressly violates time reversal and 



spin rotation symmetries, a P-odd nematic state was considered by Chubukov [0. In our 
considerations, however, the external magnetic field is zero. 

The parity operation is equivalent to interchange of r, r'. The conventional spin ne- 
matic with Qij(r,r') considered at equal times is even under P, as follows from Eq. (|I]) for 
Ajj(r, r'|£) = Qij(r, r'\t). We emphasize this distinction because it shows that the CSN state 
cannot be obtained within the old classification scheme || and is a new type of state. 

As before, we take the time derivative to write the equal-time correlator. We find 

L y (r, r') = d t Q tJ (r, r'|0) = j dr"K(r, r")e iln V lnj (r, r" , r'). (12) 

Here Vi n j(r, r", r') = (S[(r)S n (r")Sj(r')) is the three-spin correlator and an isotropic tensor 
Kij = 5{jK was taken. As mentioned above, only the P-odd part of K(r, r") contributes to 
the integral and this constraint requires special models for exchange interactions. 

The transformation properties of Vi n j under permutation of the spin indices are given by a 
particular representation of the permutation group of three objects £3 and in this particular 
case it is the I = 2 representation. The symmetry of this tensor can be seen from Eq. (|12"D: 
Vi n j = V\i n ]j, where the square brackets stand for antisymmetrization with respect to the 
indices contained within. The even-time totally symmetric three-spin correlator has been 



considered by Gor'kov [jTT|| . This symmetry is obvious from the fact that the product eu n Vi n j 
enters Eq. ([12]) • Ori the other hand, the time derivative d t Qij = Lij yields the symmetric 
traceless tensor L^. Therefore, the product eu n Vi n j should be symmetrized with respect to 
In what follows, we limit ourselves to the case of uniaxial solutions for the CSN. An 
allowed solution is 

Lij(T, r') = L(r, v'^mrij - ^8^), (13) 

where n is the nematic director (n and — n are equivalent) in spin space. The director 
specifies the plane of the spin modulation, as in the conventional spin nematic. Although 
the form of in the above equation is identical to that for the conventional even-time spin 
nematic with (r, r') = Q(r, r') {jiifij — f3j , here the physics is fundamentally different. 
The odd-time Ly is inherently connected to the three spin correlator Vi n j(r, r", r') in Eq. 
(|12|), in contrast to a two spin correlator as in the conventional spin nematic. In addition, 
the spatial parity of is opposite to the parity of the standard spin nematic Qij. 

In the exchange approximation, the coordinate and spin indices are independent and 
one can factorize V ljm : (S^r^S^i^^-^)) = T lnj $(r 1 , r 2 , r 3 ). The spin part is T lnj = 
£iin(ninj — ^Sij), and the orbital part $ is to be integrated with K{r, r') in Eq. ( [12|) to yield 
L(r, r') One finds 

Qij(r,r'\t) = Q{r,r'\t){mnj - -<%), 

d t Q(ry\t) = L(r,r'). (14) 

The lagrangian for the CSN state is identical to the lagrangian of the conventional spin 
nematics and is 

£ = ^Xii( n x 9 t n + 7 H )i( n x 9 t n + 7H)j - J kl d k ndin - U anis , (15) 

where 7 is the gyromagnetic ratio, H is the external magnetic field, is the inhomogeneous 
exchange tensor, Xij — X±(% + n i n j) is the transverse susceptibilty tensor , and U an i S is the 
relativistic anisotropy energy. Linear gradient terms are not allowed in the lagrangian. 
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We mention briefly possible ways to generate a CSN phase. To lower the symmetry of 
the system down to the CSN one can consider the quadrupolar interaction in the CSL phase 
(I = 0), which will generate the nematic (I = 2) component. In the P,T violating ground 
state this nematic state will have a chiral component. Another possibility is spontaneous 
breakdown of P and T in an already existing nematic state. 
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